Single-site entanglement at "superconductor"-insulator transition in the Hirsch model 
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We investigate the transition to the insulating state in the one-dimensional Hubbard model with bond-charge 
interaction x (Hirsch model), at half-filling and T — 0. By means of the density-matrix renormalization group 
algorithm the charge gap closure is examined by both standard finite-size scaling analysis and looking at sin- 
gularities in the derivatives of single-site entanglement. The results of the two techniques show that a quantum 
phase transition takes place at a finite Coulomb interaction u c (x) for x > 0.5. The region < u < u c turns 
out to have a superconducting nature, at least for not too large x > x c . 

PACS numbers: 71.10.Hf, 03.65.Ud, 71.10.Fd, 71.30.+h 
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I. INTRODUCTION 

Recently 1 - 2 the role of the bond-charge interaction x in in- 
ducing interesting physics in low-dimensional correlated elec- 
tron systems has been underlined. Such an interaction orig- 
inates from the Coulomb repulsion between electrons in a 
solid, and its actual value is in general smaller than the on-site 
Coulomb interaction (Hubbard u term), though not negligible 
in many cases. Moreover, the x term breaks the particle-hole 
symmetry of the model, which feature is known to be relevant 
in relation to high-T c superconductivity. In particular, in the 
one-dimensional case, the bond-charge interaction has been 
proven to play a relevant role in real materials. 3 

A complete study of the behavior of a one-dimensional 
Hubbard system in which the x term is included is still miss- 
ing, though a certain number of interesting and competing in- 
formation are known. First of all, x = corresponds to the 
standard Hubbard model, in which case at half-filling (n = 1) 
and T = a metal-insulator transition takes place at u c = 0, 
the insulating state being characterized by antiferromagnetic 
order. Moreover, bosonization studies 4 show that the critical 
value u c is independent from the actual value of the bond- 
charge interaction, at least for x <C 1. On the other hand, 
an exact solution of the x = 1 case is also available. 5 There 
the charge gap responsible for the insulating phase at n = 1 
opens at u c = 4. The nonvanishing value of u c is a clear 
indication of the failure of the bosonization approximation in 
this regime. Still, at x = 1 and for u < u c (n), the sys- 
tem enters a phase characterized by off-diagonal long-range 
order (ODLRO 6 ) and nonvanishing pairing correlation, both 
features being an indication of possible superconducting or- 
der, which is absent at x = 1 due to the high degeneracy of 
the ground state (g.s.). The possible presence of supercon- 
ducting order has been investigated also for x ^ 1 at various 
n. It turns out that this is the case, at least for low enough u 
and appropriate x and n values; 7 ' 8,9 also, a spin gap opens for 
n > l, 2 ' 8 whereas, as for the transition to the insulating phase, 
results at x ^ 1 are far from exhaustive, 910 though support- 
ing the possibility of a finite u c away from the weak-coupling 
limit. 

To resume the physical picture emerging from previous 



studies, at n = 1 and T = it is expected that the system 
undergoes a change at intermediate values of x, passing from 
a situation in which a metal-antiferromagnetic insulator tran- 
sition takes place at vanishing u, to a situation in which a pos- 
sible superconductor-insulator transition occurs at finite u. 

The purpose of the present paper is to investigate such a 
change. This is achieved by means of a composite analysis of 
the numerical results obtained using the density-matrix renor- 
malization group (DMRG) algorithm (Sec. III). On the one 
hand, the study of charge gap closure is performed by means 
of its finite-size scaling (FSS) behavior (Sec. Ill A). On the 
other hand, quite recent results 1 1 have proven that, at least at 
x = 1, the quantum phase transition (QPT) to the insulating 
state at half-filling corresponds to a divergence (in the ther- 
modynamic limit) in the partial derivative with respect to u of 
the single-site entropy of entanglement. 12 S Accordingly, we 
investigated the behavior of S also for x ^ 1 (Sec. Ill B). 
Furthermore, in Sec. Ill C both spin-gap and pair-correlation 
behaviors are analyzed within the noninsulating phase. 



II. THE HIRSCH MODEL 

The model we deal with is the bond-charge extended Hub- 
bard model first introduced by Hirsch^ It is described by the 
Hamiltonian 



H 



En 



'it "a 



(1) 



where cj CT and Ci a are fermionic creation and annihilation op- 
erators on a one-dimensional chain of length L; a =f, j is the 
spin label, a denotes its opposite, n Jcr = ct^c^ is the electron 
charge with spin a, and (i, j) stands for neighboring sites; u 
and x are the (dimensionless) on-site Coulomb repulsion and 
bond-charge interaction parameters, respectively. 

The model Hamiltonian (Q]i has su(2)®u(l) symmetry, 
since it commutes with total spin and charge operators; cor- 
respondingly, their eigenvalues 5 t z ot and N are good quantum 
numbers. Due particle-hole transformation properties, the pa- 
rameter region < x < 1 is representative of all x values. 
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In the following we shall focus on the half-filled regime. 
For x = 1 H acquires an extra symmetry, since in this case 
the so-called 7/-pairing operator 77 = J^jVj = J2j c jl c j1 
and its hermitian conjugate both commute with the Hamil- 
tonian, which also preserves the total number of doubly oc- 
cupied sites. These symmetries can be exploited to ob- 
tain the complete spectrum of the model, from which both 
zero-temperature 5 and finite-temperature 14 properties can be 
gained. In particular the g.s. phase diagram is characterized 
by two QPT's. At u = 4 a charge gap opens and the system 
enters an insulating phase, whereas for —4 < u < 4 a gap- 
less phase characterized by superconducting correlations and 
ODLRO is approached. The latter reads Yim.u_j\^ 00 (r/jrij) = 
n 2 d , with rid the average number of doubly occupied sites 
(0 < rid < 1/2). For u < —4 a further phase appears, which 
contains just empty and doubly occupied sites, and is again 
characterized by ODLRO. Noticeably, all phases containing 
singly occupied sites are highly degenerate, since the spin ori- 
entation is arbitrary. The x — 1 case has also been explored 11 
by means of single-site entanglement. In the sector with zero 
magnetization (Sfot = 0) tne latter reads 



2n d log 2 rid + (1 - 2n d ) log 2 



(2) 



In the thermodynamic limit the above QPT's appear as diver- 
gencies in d u S. In particular, the transition to the insulating 
state is signaled by an algebraic divergence of the latter quan- 
tity. S has been studied at x — as welli^ in this case no 
singularities are shown by d u S. 



III. DISCUSSION OF NUMERICAL RESULTS 



We fixed x to different values and, for each of them, we 
evaluated the charge gap (fj) for several values of u and L 
in order to be able to perform the FSS analysis and to find the 
true (that is, in the thermodynamic limit) critical value of u for 
which the gap closes. We recall that, as x grows, the system 
goes to a more symmetric and ordered phase, to the point that, 
at x = 1, the g.s. becomes highly degenerate, and for u < 4 it 
also exhibits ODLRO, which implies strong long-ranged pair 
correlations. Accordingly, in the vicinity of that region the 
DMRG numerical simulation experiences some problems; the 
program precision becomes poorer, since the truncation pro- 
cedure intrinsic to DMRG is expected to be less efficient when 
long-range correlations are present. Whereas when u > u c (x) 
— gapped phase — numerical estimations are much more ac- 
curate: this allowed us to obtain — with an acceptable error 
bar — the critical points, as reported in Fig. IHI Al which sum- 




In order to study the transition to the insulating state at arbi- 
trary x we used the DMRG method. 16 The algorithm has been 
customized 17 for the convergence of not only the g.s., but also 
of several excited-state energies. This is important when the 
g.s. in a given sector of N and S* ot is quasidegenerate and 
the procedure with only one target state could give unreliable 
results. The number of optimized DMRG states was fixed to 
(at least) m = 512, and periodic boundary conditions have 
been used. The correlation functions on the g.s. are computed 
at the end of the finite-system iterations. 



A. Charge-gap closure 

We computed the charge gap at half-filling as 

_ E(L + 2)+E(L-2)-2E(L) 
A c - 2 , (3) 

where E(N) corresponds to the g.s. energy of the system 
evaluated as the lowest-energy state in the sector S* ot = 
with N electrons. Indeed, Eq. © allows one to analyze 
the conducting properties of the system without affecting the 
magnetic ones. 



FIG. 1: (Color online) The critical values for the opening of the 
charge gap in the u, x plane, as obtained by the different meth- 
ods: exact results (triangles), FSS analysis (dots), and singularities 
in d a S, for a — u, x (crosses). M, I, and SC stand for metallic, 
insulating, and superconducting phases respectively. 

marizes our results, displaying how there are basically two 
regimes: < x < 0.5 (i) and 0.5 < x < 1 (ii). In (i), the 
transition takes place at u c — 0: a small positive u forces the 
system to become a charge insulator. Here the numerical sim- 
ulation is not so demanding, and the results of FSS analysis of 
gap closure are quite accurate, as shown in Fig. [2] 

On the other hand, in region (ii), u c becomes different from 
zero and positive, increasing smoothly to match the analyti- 
cal exact result u c (l) = 4. Correspondingly, a noninsulating 
phase is stable for positive values u < u c of the Coulomb re- 
pulsion parameter. Interestingly, the behavior of the gap clo- 
sure when entering such a phase changes; it is no longer of 
linear type, as reported in Fig. [3] 

At a; = 0.5 the model Hamiltonian acquires extra symme- 
tries, since the term interchanging singly and doubly occupied 
sites disappears. We observed that in its surrounding the nu- 
merical simulations, and in particular the FSS analysis, expe- 
rienced greater difficulties, to the point that we were not able 
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FIG. 2: (Color online) Finite-size scaling of the charge gap l[3} eval- 
uated at u = for x = 0.1, 0.2, 0.3, 0.4, and L = 10 - 62 as a 
function of 1 / L. 




FIG. 3: (Color online) Finite-size scaling of the charge gap l[3} eval- 
uated at x = 0.7, u = 2.4, 2.45, 2.5, 2.6 and L = 10 - 62 as a 
function of 1/L 2 . 



in this case to distinguish by this method whether or not the 
gap closes at finite u c (up to u c < 1). 



B. Single-site entanglement 

As already pointed out, recent results have shown how at 
least some QPT's can be spotted by measures of entangle- 
ment. In particular, the transition to the insulating phase at 
x = 1 is signaled as an algebraic divergence of <9 u iS.— We 
then perform a study of the same transition at x ^ 1 through 
d a S. Here S is a function of n c i [see Eq. an accurate 
estimate of is obtained from the correlation matrix (niUj) 
( n i — /Co- by averaging over its diagonal entries. This 
procedure compensates for the lack of a strict translational 



invariance due to the truncation approximation. Moreover, 
also finite-size effects are rather small, since S is related to 
single-site correlations and can be regarded as a thermody- 
namic quantity for L ^> 1; it turns out that its evaluation at 
L = 30 gives already a good estimate of its large-i behavior. 

The qualitative behavior of S for different x is similar, be- 
ing a decreasing monotone function of u with a flex depend- 
ing on x. Nevertheless, the behavior of its derivatives d a S 
changes drastically as x » 0.5, u > 0. In Fig. |4]we report it 
for a = u: a strong signature of the undergoing phase transi- 
tion is seen whenever x > 0.5, depicted as a sharp singularity 
at a critical point u' c (x), which should be reminiscent of the 
divergent behavior of the same quantity in the thermodynamic 
limit. Remarkably, it turns out that u c = u' c within the whole 
range of x values for which u c can be evaluated through FSS 
(see Fig. IITJ At . Moreover, for the reasons explained above, 
u' c has negligible error bar, as opposite to u c . Since S is eas- 
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FIG. 4: (Color online) d u S as a function of u, evaluated at L = 30 
for (left) x = 0.1, 0.3, 0.4, 0.5 and (right) x = 0.6, 0.7, 0.8, 0.9. 

ily worked out for arbitrary u, x, we may use the study of 
its derivatives for guessing u c , x c in those regions (x « 0.5, 
u < 1) where the study of gap closure through FSS is hard. 
Indeed, in Fig. IIII Al the points indicated just with crosses 
along the critical curve are obtained simply as singularities of 
d x S. 

Notice that the divergence in d u S disappears for x < 0.5. 
We recall that for x < 0.5 we obtained u c = 0; we then con- 
jecture u' c = = u c also at x = 0.5. We may infer that the 
different behavior in d u S is related to the phase which opens 
below the insulating phase: a smooth behavior is found when 
such phase is metallic (x < 0.5), whereas a singular behavior 
appears when the phase changes (see below). This is signaled 
by the fact that also d x S displays a singular behavior in corre- 
spondence of x c whenever u > 0. In terms of the continuum 
theory, one passes from a situation that is physically similar to 
the case x = 0, where the gap is created by a marginally rele- 
vant operator (as predicted by bosonization) and is character- 
ized by an essential singularity, to a regime where the charge 
gap A c oc (u — u c ) v (with v < oo) is generated by a truly 
relevant operator. In fact, according to the scaling analysis in 
Ref. [ill d a S should diverge with an exponent p = 2 — 2v 
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whenever v < 1, and the numerical data in the right panel of 
Fig. |4]clearly show a finite value of p and thus of v. 

As a matter of fact, let us notice that the singularities found 
in d a S are a direct consequence of the behavior of d a rid, since 
rid in the present case is the only nontrivial correlation con- 
tained in S [see Eq. (fJJ]. Hence the change which takes place 
along the critical curve is driven by rid as a function of u or x. 



C. Pairing and spin-gap 

A study of the binding energy 

E b = ^[E(L + 2)+E(L)-2E(L+l)] (4) 

confirms that for x > x c , u < u c a new phase opens, char- 
acterized by negative E b . This fact could correspond simply 
to phase separation, as well as to superconductivity, both phe- 
nomena having been observed in previous numerical studies 
by varying x, u, and n. 8,10 In order to understand the nature of 
the new phase in our case, we have studied both singlet pair 
correlations {rj\rjj) and the spin gap 



K = E{L,Sl t = l)-E{L,Sl t = Q) 



(5) 



at x = 0.6 at various L. The results support the superconduct- 
ing (SC) nature of the phase. In fact, within such phase, the 
spin gap turns out to be open up to u c (see Fig. |5). 
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FIG. 6: (Color online) Pair correlations at half chain (\i — j\ = L/2) 
vs chain length at x — 0.6 and various u values. 



we do not expect that this would differ qualitatively from that 
observed at the specific value x — 0.6, we resolved to call that 
phase SC. Notice also that in real materials it is more likely to 
have x w 0.6 than x = 1. 



IV. CONCLUSIONS 
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FIG. 5: (Color online) Spin gap vs chain length at x 
various u values. 



0.6 and 



Furthermore, in Fig. [6] it is shown that as u < u c the pair 
correlations do not vanish for \i—j\ —> oo, reaching an asymp- 
totic value of the order of n 2 ,, as in the exact case x = 1. A 
detailed study of the above phase would require an exhaus- 
tive investigation of spin gap and pair correlations for higher 
x values as well. Such analysis is beyond the purposes of the 
present paper and is expected to become more demanding as 
x becomes closer to 1. Moreover, at x — 1 the spin gap dis- 
appears. However, as far as the physics is concerned, since 



In summary, we have studied the role played by the bond- 
charge interaction in the opening of the insulating phase in 
the Hirsch model. The critical curve (u c , x c ) is accurately 
determined through the composite use of FSS analysis of 
charge-gap closure and of singularities in S; the curve qual- 
itatively agrees with results obtained within the slave-boson 
approximation. 9 The point x = 0.5 appears to play a key role. 
In particular, the bosonization description, indicating u c = 0, 
holds up to a; w 0.5, whereas the rich physics characteristic 
of the integrable (and highly degenerate) case x = 1 in fact 
qualitatively survives down to x rj 0.5. The region below u c 
changes from metallic to SC again at x « 0.5. In correspon- 
dence, the FSS analysis of gap closure shows that the depen- 
dence of A c on 1/L at the critical point passes from linear 
to nonlinear, which fact may be interpreted within conformal 
field theory 19 as a signature of vanishing charge velocity v c . 
This possibility is supported both by a direct calculation at 
u = 2.45, x = 0.70, where we found v c ~ 10~ 2 , and by 
Fig. [2] in which the line slope decreases when approaching 
x = 0.5 from below. Note that such a scenario requires that 
the central charge along the critical curve be c = 1, for both 
the charge and spin degrees of freedom. Using the formula 
in Ref. 20 for the so-called block entropy we have checked 
that this is indeed the case. In terms of the quantum-classical 
mapping 21 the 1 /L 2 behavior suggests that the dynamical ex- 
ponent is 2 for x > 0.5, which might be interpreted as the 
system approaching a parabolic point in the band. 

The procedure of using the single-site entanglement as a 
witness of ongoing QPT's is quite convenient if considered 
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from a time-consuming perspective, since it does not require 
a thorough FSS analysis (at least in this context). More- 
over, once it is found that a particular transition (like the SC- 
insulator transition) is captured by such procedure, this could 
be used beyond the one-dimensional case, by simple imple- 
mentation of exact-diagonalization techniques on small two- 
dimensional clusters. 

Even within the one-dimensional case, moving away from 
half-filling should provide other interesting physics, since 
particle-hole invariance for x ^ 1 is lost. It is expected that 
within the hole-doped region higher values of u c for the SC 
phase could be achieved; 8 though such a QPT is of a different 



nature at n ^ 1, it should be captured by appropriate measures 
of multipartite entanglement^ 
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